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Abstract 

We discuss stability of dynamically collapsing gas spheres. We use a similarity solution for a dynamically 
collapsing sphere as the unperturbed state. In the similarity solution the gas pressure is approximated by 
a polytrope of P = Kp'^. We examine three types of perturbations: bar {£ = 2) mode, spin-up mode, and 
Ori-Piran mode. When 7 < 1.097, it is unstable against bar-mode. It is unstable against spin- up mode for 
any 7. When 7 < 0.961, the similarity solution is unstable against Ori-Piran mode. The unstable mode 
grows in proportion to \t — to| while the central density increases in proportion to pc (x {t — to)^^ in 
the similarity solution. The growth rate, a is obtained numerically as a function of 7 for bar mode and 
Ori-Piran mode. The growth rate of the bar mode is larger for a smaller 7. The spin-up mode has the 
growth rate of a = 1/3 for any 7. 
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1. Introduction 

As well as the majority of nearby pre-main sequence 
stars, most young stars have their companions (see, e.g., 
Mathieu 1994). Since even young protostars have their 
companions, fragmentation is more plausible for forma- 
tion of multiple star systems than capture and other pro- 
cesses. 

Fragmentation during the star formation has been paid 
much attention and a number of numerical simulations 
have been performed for clarifying it. Nonetheless the 
results are not conclusive. Some early simulations have 
only limited spatial resolution and their results on frag- 
mentation are unreliable (Truelove et al. 1997). A gas 
cloud fragments artificially in numerical simulations of 
low resolution. If we exclude the artificial fragmentation, 
fragmentation proceeds but only slowly in numerical sim- 
ulations. Fragmentation competes with the total collapse 
of a cloud. A cloud collapses and the density increases 
before it fragments. Only when an initial model was very 
much elongated, the simulation showed fragmentation of 
a cloud at a relatively low density. This initial model is, 
however, not applicable to a roundish dense core contain- 
ing a binary. 

Fragmentation has been studied not only by numerical 
simulations but also by linear stability analyses. While 
numerical simulations are advantageous at handling non- 
linearity, linear stability analyses are helpful for elucidat- 
ing underlying physics. In this sense these two technolo- 



gies are complementary. Most linear stability analyese 
thus far, however, deal with the stability of an equilib- 
rium cloud against small perturbation. The competition 
with the total collapse could not be taken account into 
them. Silk, Suto (1988) tried to compute the stability 
of a similarity solution for the collapse of an isothermal 
gas cloud but could not obtain a self-consistent solution. 
Only recently Hanawa, Matsumoto (1999) succeeded in 
analyzing the stability of a collapsing isothermal cloud 
against non-spherical perturbations. They found that a 
collapsing isothermal cloud is unstable against bar mode 
{£ — 2) perturbation. The bar mode grows in proportion 



to p, 



0.354 



where pc denotes the central density. When 



it grows, a collapsing gas cloud may change its shape 
from sphere to filament. This result is consistent with 
recent three-dimensional simulations by Truelove et al. 
(1998) and Matsumoto, Hanawa (1999) in which a fil- 
ament forms in a collapsing cloud. The former showed 
that the filament fragments to form a binary. Thus the 
bar mode is interesting for understanding fragmentation 
during the collapse. 

In this paper we extend the linear stability analy- 
ses of Hanawa, Matsumoto (1999) for collapsing non- 
isothermal spheres. For simplicity we employ the poly- 
tropic relation, P = Kp'*, for the model cloud. The 
polytrope model can describe temperature change dur- 
ing the collapse in the most simple form. The polytrope 
model of 7 ~ 4/3 can be appHed also to a collapsing 
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iron core resulting Type II supernova. In section 2 we 
present the similarity solution for gravitationally collaps- 
ing sphere of a polytropic gas cloud. This section is essen- 
tially the review of Yahil (1983) and Suto, Silk (1991). 
We summarize their results to investigate the stability. 
In section 3 we examine the stability of the similarity 
solution against bar mode. The growth rate of the per- 
turbation is larger when 7 is smaller. The bar mode is 
degenerate and its growth rate is independent of the az- 
imuthal wavenumber, m. The shape of the core suffering 
the bar mode depends on m and the sign of the pertur- 
bation. We show it visually. In section 4 we show that a 
spherically collapsing gas cloud is unstable against spin- 
up. This result is an extention of the stability analysis of 
Hanawa, Nakayama (1997). In section 6 we show that the 
similarity solution is unstable against a spherical pertur- 
bation when 7 < 0.961. The analysis shown in section 5 
is based on Ori, Piran (1988). A short summary is given 
in section 6. 

2. Similarity Solution 

For simplicity we consider gas of which equation of 
state is expressed by polytrope, 



P = Kp-f , 



(1) 



where P and p denote the pressure and density, respec- 
tively. The hydrodynamical equations are then expressed 
as 



| + v.(,.) = o 



(2) 



and 



dt 



(pv) + VP + V ■ {pv ®v) + = , (3) 



where v and $ denote the velocity and gravitational po- 
tential, respectively. The gravitational potential is re- 
lated with the density distribution by the Poisson equa- 
tion, 



A<I> AnGp , 



(4) 



where G denotes the gravitational constant. 

For later convenience, we introduce the zooming coor- 
dinates of Bouquet et al. (1985) to solve equations (|^) 
through (|). The zooming coordinates, {$,, r), are related 
with the ordinary coordinates, (r, t), by 



Co \t - to\ 
- In |1 - t/to\ 



(5) 



where co denotes a standard sound speed and is a func- 
tion of time t. The symbol, to, denotes an epoch at the 
instant of the protostar formation. The density in the 



zooming coordinates, g, is related with that in the ordi- 
nary coordinates, p, by 



g{x, t) = AttG p(t - to)'' 



(6) 



We define the standard sound speed, Cq, so that it denotes 
the sound speed at a given t when g = 1. Thus it is 
expressed as 



Co 



(7) 



The pressure in the zooming coordinates, p, is related 
with the that in the ordinary coordinates, P, by 



P 



AttG 



Pit ~ to)2 



(8) 



Substituting equations (^) and (g|) into equation (|l|), we 
obtain the polytrope relation in the zooming coordinates. 



P 



1 



(9) 



The velocity in the zooming coordinates, it, is defined as 
V ^ r 



w = - + (2 - 7) -t; rr 

Co Co \t - to 



(10) 



This velocity denotes that with respect to the zooming 
coordinates, and includes the apparent motion, the last 
term in equation (p^. The gravitational potential in 
the zooming coordinates, 0, is related with that in the 
ordinary coordinates, $ by 



(11) 



In the zooming coordinates, the hydrodyanmical equa- 
tions are expressed as 



dg_ 
dr 



+ V • {gu) = (4 - 37) ^. , 



(12) 



— {gu) + V' • [gu ®u) + V'p + (?V> 

OT 

= (2 - 7) (7 - 1) gi + {7 - 57) gil?,) 



and 



(14) 



for t < to- The symbols, V and A', denote the gradient 
and Laplacian in the ^-space, respectively. 

Assuming stationary in the zooming coordinates 
{d/dr = 0) and the spherical symmetry {d/dO = 
d/d^p ~ 0), we seek a similarity solution. Under these 
assumptions equations (|l^), ([l3|), and (|l^) reduce to 



dur 



Ur dg 

7 W 



(4 - 37) 



2Ur 



(15) 



No. 



dUr 
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= (2 - 7) (7 - 1)C 



1 / dp\ dg d4> 

~g [d~g) + a? 

+ (3 - 27) Ur , (16) 



and 



e 



gUr 



4 - 37 ' 



(17) 



where ^ — |^|. After some algebra we can rewrite equa- 
tions (|l|) and (^ into 



,7-1 



gUr 

4-37 



+ (2 - 7)(7 - l)e + (7 - l)Ur + 



2u. 



21 



(18) 



and 



(e 



dUr 

(2 - 7) (7 - 1)?". 
+ (4 - 37) f?^-^ 



4-37 

(3 - 27)7/^2 

n7-l 



2Ur 



(19) 



Equations (|T^) and (|T^) are singular at the sonic point, 
Ur"^ = g'^ ~ .We obtain the the similarity solution by in- 
tegrating equations ( p^ ) and (|l^) with the Runge-Kutta 
method. In the numerical integration we used the auxil- 
iary variable of Whitworth, Summers (1985), s, defined 

by 



(20) 



Using equation (^), we rewrite equations (|8|) and ( |ig| ) 
into 



dg 
ds 



Q 
(7 



gur 



37 



1) Ur 



-I- (2 - 7)(7 - 1)C 



2Ur 



(21) 



and 



ds 



gur^ 



37 



(2 - 7) (7 - l)e"r. 



(3 - 27) + (4 - 37) g 

T 



,7-1 



^7- 1 



(22) 



respectively. 

Similarity solutions exist for 7 < 4/3. Figure 1 shows 
the similarity solutions for 7 — 0.9, 1.0, and 1.1. The 
solid curves denote g while the dashed curves denote the 
infall velocity, — m -|- (2 — 7) ^. These solutions are the 
same as those obtained by Yahil (1983) and Suto, Silk 
(1988). 




Fig. 1.. The similarity solutions are shown for 7 = 0.9, 1.0, 
and 1.1. The solid curves denote the density, q. The dashed 
curves denote the infall velocity, — « -|- (2 — 7) ^. The infall 
velocity does not include the apparent motion in the zooming 
coordinates and is positive for inward flow. 



3. Bar Mode 

In this section we consider a non-spherical perturba- 
tion around the similarity solution. The density per- 
turbation is assumed to be proportional to the spherical 
harmonics, Yp{9, ip). Then the density and velocity are 
expressed as 



g = go + Sg{Oe--Yr{e,ip) , 



ue = (5ue(0 



d 



1 86 



(23) 
(24) 

(25) 



d 



(i + 1) sin 9 dip 



Yr{0,p), (26) 
(27) 



where the symbols with suffix, 0, denote the values in 
the similarity solution and the symbols with the symbol, 
S, denote the perturbations. Substituting equations ( [2^ ) 
throughout ( p7[ ) into equations (p^), (|l3[), and (p^), we 
obtain the perturbation equations, 

1 d 



goSue 







(28) 



(cr + 27 - 3) Sur + 



9 , J. , d 

— {UroSUr) + ^ 
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[Vol. 




+ ST = 0, 



(29) 



(cr + 27 - 3) Sue 
e + 1 Sg 



UrO d 



ST, 



2-7 



= 



and 



d 2ST 



£{£ 



■ S(f) + Sg 



(30) 



(31) 



(32) 



These perturbation equations have singularities at the 
origin (^ = 0), the sonic point [(uro)^ ~ ^ and the 
infinity (^ = +cx)). These perturbation equations reduce 
to those of Hanawa, Matsumoto (1999) when 7 = 1. 
We solve these perturbation equations as an eigenvalue 
problem for the growth rate, tr, according to them. The 
details of the numerical procedures are given in Hanawa, 
Matsumoto (1999). The growth rate is obtained as a 
function of 7 and £. The growth is independent of m 
since the unperturbed state is spherically symmetric (see, 
e.g., Hanawa, Matsumoto 1999). 

Figure 2 shows the growth rate, cr, for £ = 2 mode as 
a function of 7. The open squares denote the numeri- 
cally obtained growth rates and the curve labeled "Bar" 
denotes a smooth fit to them. The growth rate is larger 
for smaller 7. A collapsing gas sphere is unstable only 
when 7 < 1.097. We cannot find a damped mode since 
appropriate boundary condition are not given at ^ = 00 
for it (see appendix 2). 

We searched for £ — "i mode but could find none. Al- 
though £ = 1 mode exists, it denotes only the misfit 
to the center of gravity and is not relevant to real in- 
stability (Hanawa, Matsumoto 1999). A higher £ mode 
is also unlikely to exist. Thus only the bar mode gives 
a non-spherical density perturbation growing during the 
collapse. Since the bar mode is stabilized for 7 > 1.097, 
a similarity solution for 7 > 1.097 is stable against any 
non-spherical density perturbation. This result can an- 
swer to the question posed by Goldreich, Lai, Sahrling 
(1997). They asked whether a dynamically collapsing 
iron resuting type II super nova is unstable against a 
large scale non-spherical perturbation. Since the dynam- 
ically collapsing iron core is well approximated by a sim- 
ilarity solution for 7 ^ 4/3, it is stable against any non- 
spherical density perturbation. 

Figure 3 illustrates the variety of the bar mode. Each 
panel shows a dynamically collapsing isothermal gas 
cloud suffering a different bar mode perturbation. It 



t> 0.2 
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Fig. 2.. The growth rate, a, is shown as a function of 7. The open 
squares denote numerically obtained data points. 



shows the density distribution perspectively by isoden- 
sity surfaces. Panel (a) denotes a bar mode of to = 0. 
In panel (a) the collapsing gas cloud is elongated in the 
z-direction. Panel (b) denotes the same bar mode of 
m = but having the opposite sign. In panel (b) the 
collapsing gas cloud is oblate and compressed in the z- 
direction. Panels (c) and (d) denote the bar modes of 
TO = 1 and 2, respectively. When the modes grow, the 
collapsing gas cloud becomes triaxial. 

4. Spin-up Mode 

In this section we consider the velocity perturbation 
expressed as 




^ saiO dip 



(33) 



^ de 

Substituting equation ( ^3[ ) into equation ( p2[ ) we obtain 
Sg = 0, (34) 

for this mode. Similarly we obtain 

S(j) = 0, (35) 

by substituting equation (Sj) into equation (p^. Substi- 
tuting equations (p3h, (pi), and (p9) into equation (O), 
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we obtain 



Fig. 3.. A collapsing isothermal cloud core suffering the bar mode. 
The density distribution is shown by the isodensity surfaces. 
The velocity in the x — z plane is shown by arrows. Each 
panel shows a different bar mode. Panel (a) denotes the bar 
mode of m = 0. Panel (b) denotes the same bar mode of 
m = but having the opposite sign. Panel (c) denotes the 
bar mode of m = 1 while panel (d) does that of m = 2. 
In panels (a) and (b) the perturbation is added so that the 
radial velocity is Vr oc [—2/3 ± 0.10(3cos^S — 1)] r near 
the center, respectiveley. In panel (c) the radial velocity is 
Vr OC (—2/3 + 0.05 sin 29 cos (p) r near the center. In panel (d) 
it is Vr oc (—2/3 + 0.10 sin^ cos2i^) r near the center. 



((T-7 + 57) goAi^m + 



1 d 



Substituting equation (nSl) into equation (|3 



27) 



Substituting equation (A.2) into equation ( 



{S,'^ QQUroAl^m) = .(36) 

we obtain 
. (37) 
\\ we obtain 



d . , 




i = o 



.(38) 
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The growth rate, cr, should be smaller than 1/3 since 
the velocity perturbation is regular at the origin only 
for d\nAg^rn/dln£^ > 2 [see equation (^]. When 
dlnAf m/dln^ = 2, the angular velocity is nearly con- 
stant around the origin. The growth rate of the spin-up 
mode, CT = 1/3 is also shown in Figure 2 for comparison 
with that of the bar mode. This spin-up mode is essen- 
tially the same as the spin-up mode shown in Appendix 
of Hanawa, Nakayama (1997). 

5. Ori-Piran Mode 

In this section we consider a perturbation emanating 
from the sonic point. According to Ori, Piran (1988) we 
analyze a spherical perturbation which vanishes inside 
the sonic point. Substituting £ = into equations (I2J 
through (|^ we obtain 



where equations (p2|) and (43) are substituted. Equation 



„ duro - . a In go 

2 + (7 - 1) UrO ■ 



7-57 



.(46) 



When 7 = 1, equation ( [iq ) is equivalent to equation 
(15). To elucidate the growth (or damping) of this mode 
we rewrite equation (Hfl) into 



-2 — (Uro 



Cs) + 



7-57 



(47) 



d_ 



Sg - 



ieiSguro + QoSur)] = (4-37)5^,(39) 



— dUr 

OT 



6Ur 



dUrO 

6T = {3 - 27) dUr , 



Ug] 
(40) 



and 



6r 



1 

e 



6giC) e dC 



We evaluate equations (|39|), (| 
point. Then we obtain 

d 



(41) 

and (|4^) at the sonic 



UrO-^Sg = 



S(j) = 0, 



(42) 

(43) 
(44) 



at the sonic point. Equations ( [42|) and (|4^) are equivalent 
and indicate that the sound wave traveling outwards with 
the phase speed Ur + Cg vanishes. In other words the 
perturbation emanating from the sonic point is the other 
sound wave of which phase velocity is zero at the sonic 
pint, Ur ~ Cs = 0. Taking the linear combination of the 
.^-derivatives of equations (|3^) and (^0|) and evaluating it 
at the sonic point, we obtain 



d_ 



^ A 



^ A 



-2 



dUrQ 



-I- (7 - 1) UrO 



d In gio 



7-57 



' ^A ^A 
"rO ^Og - go —OUr 



(45) 



The first term in the right hand side of equation (^) de- 
notes the dilution of the wave due to the spatial variation 
of the phase velocity. The phase speed vanishes at the 
sonic point and increases with the radial distance. Thus 
the wave dilutes and the amplitude decreases. The sec- 
ond term in equation (^) denotes the self-reproduction. 
As shown in the right hand side of equations ( p^ ) and 
(|l3|), the density and momentum density reproduce and 
amplify themselves when measured in the zooming coor- 
dinates. When the self-reproduction overcomes the di- 
lution, the collapsing gas sphere is unstable against the 
Ori-Piran mode. 

Nonlinear growth of the Ori-Piran mode has not been 
studied yet. This is the first report that a similarity solu- 
tion of collapsing gas sphere can be unstable against the 
Ori-Piran mode. As far as we know, Ogino, Tomisaka, 
Nakamura (1999) is the first report on coUaps of a poly- 
tropic gas sphere with 7 < 1. They reported a little 
on their numerical simulations of 7 < 1 and mentioned 
nothing on the Ori-Piran instability. Since the similarity 
solution is unstable, their numerical solution may not be 
well approximated by it. 

6. Summary 

As shown in the previous sections a gravitationally col- 
lapsing polytropic gas sphere can suffer from three types 
of instability. When 7 < 4/3, it is unstable against spin- 
up mode shown in section 5. The spin-up mode grows in 
proportion to \t — tQ\~^^^ and accordingly in proportion 
to p^/^. When 7 < 1.097, it is unstable also against the 
bar mode instability. The growth rate of the bar mode 
is large for smaller 7. When 7 < 1.006, it is larger than 
that of the spin-up mode. When 7 < 0.961, a gravita- 
tionally collapsing polytropic gas sphere is unstable also 
against the Ori-Piran mode. 

It is particularly interesting that the bar mode has a 
larger growth rate when 7 is smaller. This means that the 
bar mode grows faster when the sound speed decreases 
with the increase in the density. This result may be valid 
also for a gravitationally collapsing gas sphere in which 
turbulence has an effective pressure. Since the effective 
sound speed is lower in a denser part of molecular cloud, 



No. 
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a collapsing cloud core may be more unstable against 
the bar mode than a purely isothermal gas sphere. This 
enhancement in the bar mode instability may be relevant 
to fragmentation of a collapsing cloud core and formation 
of multiple stars during the collapse. 
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Appendix 1. Similarity Solution for Collapse of 
Polytropic Gas Sphere 

In this appendix we summarize the main characteris- 
tics of the similarity solution for collapse of polytropic 
gas sphere. See Yahil (1983) and Suto, Silk (1988) for 
the derivation and more details. 

When we expand the similarity solution around ^ = 
by the Taylor series, they are expressed as 



eo(0 = £^0(0) - 



MO)]' 



67 



eo(0) 



(Al) 



and 



(2-7)-3 



bo(0)] 



1-7 



157 



0(e^).(A2) 



In the region of ^ 1 the similarity solution has the 
asymptotic form of 



and 



[Ur - (2 - 7) e] CX e 



(1-7)7(2-7) 



(A3) 



(A4) 



Appendix 2. Asymptotic Behavior of Bar Mode 
Perturbation 

We can derive the asymptotic forms for the bar mode 
perturbation from the requirement that the perturbation 
is regular at ^ = 0. The density and velocity perturba- 
tions around ^ = are expressed as 



Sg 
6ur 
Sue 



pee-' 
Pie + i)^-' 

- { a + [cr + 27 



(A5) 
(A6) 
(A7) 



+ 



e{4 ~ 37) 



(A8) 



where a and f3 are free parameters. When we derive 
the above Taylor series expansion, we use equation ( A2 ) . 
The derivation is essentially the same as that shown in 
Hanawa, Matsumoto (1999). 

As a boundary condition we assume that the relative 
density perturbation, 6g/ go, is vanishingly small at infin- 
ity, ^ — 00. After some algebra we obtain the asymptotic 
relations, 



Qo 



SUr OC ^-(- + 7-1)7(2-7) 

Sue oc ^-(- + 7-1)7(2-7) 

oc ^-(--27 + 2)7(2-7) ^ 



and 



(a - 27 + 2){a -37 + 4) 



(2 - 7)= 



r 6g 



(A9) 

(AlO) 
(All) 
(A12) 

1 ~^ 

£(1+1) 

(A13) 



When we derive the above relations, we use equations 



^) and (A4). See also Hanawa, Matsumoto (1999) for 
the derivation. 
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